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ABSTRACT 

The cosniological star formation rate in the combined Chandra Deep Fields North and South is derived 
from our X-Ray Luminosity Function for Galaxies in these Deep Fields. Mild evolution is seen up to 
rcdshift order unity with SFR ~ (1 + z) 2 ' 1 ■ This is the first directly observed normal star-forming 
galaxy X-ray luminosity function (XLF) at cosmologically interesting redshifts (z>0). This provides 
the most direct measure yet of the X-ray derived cosmic star-formation history of the Universe. 
We make use of Bayesian statistical methods to classify the galaxies and the two types of AGN, finding 
the most useful discriminators to be the X-ray luminosity, X-ray hardness ratio, and X-ray to optical 
flux ratio. There is some residual AGN contamination in the sample at the bright end of the lu- 
minosity function. Incompleteness slightly flattens the XLF at the faint end of the luminosity function. 

The XLF has a lognormal distribution and agrees well with the radio and infrared luminosity 
functions. However, the XLF does not agree with the Schechter luminosity function for the HoLF 
indicating that, as discussed in the text, additional and different physical processes may be involved 
in the establishment of the lognormal form of the XLF. 

The agreement of our star formation history points with the other star formation determinations in 
different wavebands (IR, Radio, Ho) gives an interesting constraint on the IMF. The X-ray emission in 
the Chandra band is most likely due to binary stars although X-ray emission from non-stellar sources 
(e.g., intermediate-mass black holes and/or low-luminosity AGN) remain a possibility. Under the 
assumption it is binary stars the overall consistency and correlations between single star effects and 
binary star effects indicate that not only is the one parameter IMF(M) constant but also the bivariate 
IMF(Mi, M 2 ) must be constant at least at the high mass end. Another way to put this, quite simply, 
is that X-ray observations may be measuring directly the binary star formation history of the Universe. 

X-ray studies will continue to be useful for probing the star formation history of the universe by 
avoiding problems of obscuration. Star formation may therefore be measured in more detail by deep 
surveys with future x-ray missions. 

Subject headings: galaxies, cosmology, star formation, surveys, x-rays 
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There have been many recent studies of star for- 
mation in galaxies and of the star formation his- 
tory of the universe derived from da t a in the ra- 
dio, IR, and optical (iLillv et alJ 119961 iMadau et alJ 



1998HRowan- Robinson et alJll997t|Haarsma et alJl200ft 
" Cole et al.H200lHBaldrv et alJl2002HTeplitz et alJl2003|) . 

In the range from the present epoch to redshifts of 
order un ity, recent cri t ical c o mpilations and discus - 



_ ^ , iHnn kius et al. ( 2003) , 
Hoed (2004) show that the re- 



signs of iSulliyan et alJ . 
Sulli van et all l|2004|) and 
suits from the multi-waveband data have a dispersion 
of 1-2 orders of magnitude in the comoving cosmic star 
formation density . As noted by these authors, there 
are important physical corrections that need to be made 
to go from the observations in a particular band to the 
cosmic star formation rate which include physical under- 
standing of the dust extinction, the Initial Mass function 
(IMF), and stellar population models. Reasonable in- 
terpretations of the current observations in this redshift 
range have been presented that argue, on the one hand, 
that there is a ste ep dec l ine in the star formation rate to 
the present epoch|Hogg (2004) or, on the other, that the 
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cosmic star formation density has a shallow decline in 
the sa me redshift range IWilson. Cowie. Barger fc Burke 
(2002). Therefore it is important to utilize all wave- 
bands to study this phenomenon from different aspects 
and with different selection effects. The X-ray band has 
now just opened up to detailed studies of star formation 
in galaxies at cosmological distances. 

Hitherto, even deep X-ray surveys studied only the 
cosmological populations of evolving active galaxies and 
quasars. X-ray studies of individual nearby galaxies were 
performed and the underlyin g hot gas and s tellar x-ray 
source components analyzed (lFabbianolll989|) . However, 
the deep 1-2 Megasecond surveys in the Chandra Deep 
Field South and North, respectively, now show a major 
cosmological population (in the range from present day 
to redshifts of order unity) of X-ray emitting normal star 
forming galaxies at faint flux levels. Stacking analysis al- 
lows one to extend th e x-ray properties studied t o even 
larger redshift ranges (JHornschemeier et al J 120 02'). 

Some notable results on star-forming galaxies have 
already been derived in deep X-ray survey work. In 
the Chandra Deep Field-North (hereafter CDF-N) it 
has been found that the faint 15^m population, com- 
posed pri marily of lumin o us inf rared starburst galax- 
ies (e.g. iCharv fe Elbazl l)2001|) ) are associated with 
X-ray-detected emission-line galaxies (Alexander et al. 
2002). Chandra and XMM-Newton stacking analyses of 
relatively quiescent (non-starburst) galaxies have con- 
strained the evolution of X-ray emission with respect to 
optical emission to at most a factor of 2-3 (e.g., Horn- 
schemeier et al. 2002; Georgakakis et al. 2003). Stud- 
ies of the quiescent population of galaxies has also pro- 
vided some initial constraints on the evolution of star- 
formation in the Universe ( SFR ~ (1 + z) k for z « 1 
where k < 3; Georgakakis et al. 2003). 

We now discuss some relevant details of the X-ray 
observations. We derive X-ray Luminosity Functions 
and the corresponding cosmological star formation his- 
tory using X-ray data from the 2 Ms and 1 Ms ex- 
posures of the CDF-N and Chandra Deep Field South 
(hereafter CDF-S). This analysis requires the extreme 
depth of the CDF surveys as non-active galaxies arise 
in appreciable numbers at extremely faint X-ray fluxes; 
the fraction of X-ray sources that are normal and star- 
forming galaxies rises sharply below 0.5-2 keV fluxes 
of » 1 x 10~ 15 erg cm " 2 s _1 (e.g., Figure 6 of 
IHornschemeier et al.l2003l ; see also the fluctuation anal- 
ysis results of Miyaji & Griffiths 2002). 

The luminosity functions for galaxies in the radio, op- 
tical and in f rared have been studied extensively (c.f. 
iTresse et alJ (|2002|) for a recent discussion). Prior to 
this work, the X-ray luminosity function for normal 
galaxies was estimated from the optical galaxy lumi- 
nosity function bv lGeorgantopoulus. Basilakos fe Plionisl 
(1999). Schmidt, Boiler & Voges (1996) reported on a 
galaxy XLF derived in the local (v < 500 km s _1 ) Uni- 
verse as well, however their sample does include some 
fairly active galaxies (e.g., Cen A) and does not appear 
to be an entirely clean normal/star- forming galaxy lu- 
minosity function. Our work seeks to construct a rela- 
tively clean normal star-forming galaxy luminosity func- 
tion, and possesses the great advantage of doing this at 
cosmologically interesting distances (our two luminosity 
functions have median redshifts z ~ 0.3 and z « 0.8). 



We may thus study the evolution of X-ray emission from 
galaxies with respect to cosmologically varying quantities 
such as the global star-formation density of the Universe. 
Our galaxy XLFs have their own set of selection effects 
that we discuss later but they can add to the overall pic- 
ture of the physics of galaxy evolution. For example, it 
is obvious that corrections for obscuration will not be 
as important here as in the analysis of Ha Luminosity 
functions discussed in detail below. 

We expect correlations between the various classic cos- 
mological indicators of star formation (Ha, IR, radio) 
and our X-ray studies. In general, the fluxes we measure 
in the respective wavebands associated with the differ- 
ent methods are all ultimately connected with the evo- 
lution and death and transfiguration of massive stars. 
For example, the radiative luminosity of OB stars, the 
mechanical energy input from supernovae and the ac- 
cretion power of high mass X-ray binaries (HMXBs) are 
all ultimately derived from massive stars. It is therefore 
not surprising that tight correlations have been found 
empirically between the X-ray flux, the radio flux and 
the IR continuum. The mo st recent careful analysis 
in the local Universe is from IRanalli. Comastri fc Settl 
(2003) which is discussed in detail below. Interestingly, 
the correlations appear to hold fo r at least the radi o 
and X-ray bands at high redshift IBauer et al.l {2002); 
IGrimm. Gilfanov fc Sunvaevl l|2003j) . We use these latest 
empirically derived correlations betwen radio, IR, and X- 
ray star formation indicators as the central relations that 
allow us to infer star formation rates from the X-ray data 
on galaxies in the Chandra Deep Fields North and South. 
In addition, it is useful to also consider theoretical stud- 
ies on the cosmic X-ray evolution of galaxies whi c h hav e 
been presented b y ICayaliere. Giacconi fc Mencil (2000), 
IPtak et all J200l and lGhosh fc White! (I200lfr 

Observations of detailed X-ray stellar populations 
in individual nearby galax i es have been c arried out 
bv iZezas fc Fabbianol l |2002|)lKilgard et al. I l)2002|) and 
iColbert et alJ l)2003|) . Using archival Chandra data on 
32 galaxies in the nearby Universe, Colbert et al. (2003) 
have established that the X-ray emission from accreting 
binaries in galaxies is correlated with the current level 
of star formation in their host galaxies. This relation 
appears to be both macroscopic in that the entire X- 
ray point source luminosity of galaxies scales with star- 
formation rate and microscopic in that the slope of the 
X-ray binary luminosity function within galaxies corre- 
lates with star formation. 

In our studies, much care has gone into accurately 
selecting the population of normal star forming galax- 
ies. Contamination by AGNs is a potentially serious 
problem. We have a multi-parameter probability space 
for the selection of the normal star forming galaxies 
and we have therefore used Bayesian methods to clas- 
sify the galaxies and the two types of AGN. The most 
useful discriminators are the X-ray luminosity and X- 
ray hardness ratio as well as the optical to X-ray ra- 
tios. AGN contamination is most serious at the bright 
end of the luminosity function. Our analysis is a step 
towards a fully modern statistical Monte Carlo Markov- 
Chain study but we need to have a better idea of the 
overall structure of the probability space be f ore w e un- 
dertake this (c.f. iHobson. Bridle and Lahavl l)2002[) and 
iHobson fc McLachlanl l|2003D L Our analysis assumes, as 
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a first step, that the priors are Gaussian and we show 
that this is a reasonable and useful assumption. 

These studies of the cosmological evolution of normal 
star forming galaxies in the x-ray bands are the first of 
many detailed studies that can be undertaken with future 
missions such as XEUS. Here we have made a start on 
the basic luminosity functions and cosmic star formation 
histories. 

The organization of this Paper is as follows: in section 
2 we describe the data acquisition, the selection of the 
galaxy sample and the processing of the data. Section 3 
presents the derivation of the galaxy XLF. In Section 4, 
we carefully compare the galaxy XLF derived here with 
the I R LF using the recent empiricall y derived analysis 
from lRanalli. Comastri fc Settil l)2003j) (hereafter RCS). 
A similar procedure is undertaken for comparison with 
the HoLF. We also give the derived cosmic star forma- 
tion rates from the XLF. In section 5 we discuss the im- 
plications of our results and the potential consequences 
for future missions. In section 6 we give our four prin- 
cipal conclusions. Throughout this paper, if not explic- 
itly stated otherwise, we use H = 70 km s" 1 Mpc -1 , 
fl m = 0.3 and fi A = 0.7. 

2. DATA ACQUISITION, SELECTION OF GALAXIES, AND 
DATA REDUCTION 

2.1. Data Acquisition 

CDF-S 

From October 1999 to Dec 2000, 11 individual expo- 
sures of the CDF-S were performed with the Chandra 
ACIS instrument, resulting in a 1 Ms exposure. The sen- 
sitivity of this deep exposure has reached flux limits of 
5.5 x 10~ 17 erg cm -2 s _1 in the soft band (0.5-2 keV) and 
4.5 x 10 -16 erg cm~ 2 s _1 in the hard band (2-10 keV). At 
these flux levels, > 80% of the cosmic X-ray background 
in both b ands is resolved. A total of 346 sources has been 
detected ijRosati et alJ 12003 IGiacconi et alJ l2002|) : the 
X-ray data reduction, so urce detection, a nd X- ray source 
catalog can be found in IGiacconi et alJ ([2002) . The op- 
tical spectroscopic follow-up observations were obtained 
using FORS1/FORS2 on the VLT for the possible optical 
counterparts of 238 X-ray sources, yielding spectroscopic 
redshifts for 141 X-ray sources. The optical spectr a 
and the redshifts are presented in ISzokolv et al.l {2004). 
Zheng et al., in preparation, have used the ten near-UV, 
optical, and near-infrared bands to estimate photomet- 
ric redshifts for 342 (99%) of the CDF-S sources, making 
detailed comparison with the spectroscopic redshifts. In 
some cases, the optical spectroscopic redshifts of Szokoly 
et al. were not highly confident and Zheng et al. was 
able to place a better redshift constraint using the multi- 
wavelength SED of the sources. We therefore used those 
redshifts of Zheng et al. which we consider to supercede 
those of Szokoly et al. 

CDF-N 

From November 1999 to February 2002, 20 individ- 
ual exposures of the CDF-N were performed by the 
Chandra ACIS instrument, resulting in a 2 Ms expo- 
sure. In the central parts of the CDF-N the sensitiv- 
ity reaches 2.5 x 10~ 17 erg cm -2 s _1 in the 0.5-2 keV 
band and 1.4 x 10~ 16 erg cm~ 2 s _1 in the 2-8 keV 
band, the 1 Ms survey; 503 sources were detected in 



the 2 Ms data, and the X-ray data reduction, source 
detection, and X-ray source catalog can be found in 
I Alexander et alJ (|2003). The optical spectroscopic and 
photometric follow- up observations were obt ained using 
Keck and Subaru (JBarger et al J 120021 12003). Spectro- 
scopic redshifts have been obtained for 284 of CDF-N 
X-ray sources; using the multiwavelength photometric 
data, photometr ic redshifts w e re ob tained for an addi- 
tional 78 sources iBarger et alJ J2003). We use these 362 
redshifts in our study. Note that although the overall 
completeness is roughly 71%, including the photomet- 
ric redshifts, the spectroscopic completeness alone for R 
< 24 is 87%. The optical and near- infrared photom- 
etry, spectrosco pic redshifts, and p hotometric redshifts 
are presented in IBarger et al.l ((2003). 

2.2. Galaxy Selection 

In our study it is essential to classify galaxies cor- 
rectly. We have employed two different galaxy classifi- 
cation techniques. The first method uses standard op- 
tical classification for sources with high-quality optical 
spectra. The second approach utilizes a multi-parameter 
Bayesian analysis. 

We have chosen to construct XLFs in the 0.5-2 keV 
band since for moderately obscured AGN with signif- 
icant star formation present, the s oft band is dom- 
inated by star forming proc e sses (IPta k et alJ 119991 
i Levenson. Weaver. &: Heckmanl 120011 iTerashima et alJ 
120021) . 

Optical Spectroscopic Classification 

Very few galaxies have optical spectra with good 
enough signal-to-noise to allow classifications based on 
emission-line flux ratios. Therefore, it is essential that 
we develop a classification scheme which is not dependent 
(or overly dependent) upon optical spectroscopic classifi- 
cation. However, there is important optical spectroscopic 
information contained within the spectra. For the CDF- 
S spectra with sufficient si gnal-to-noise, we have used 
the di agnostic diagrams of Rol a. Te rlevich k. Terlevichl 
(1997) in a manner similiar to iKewlev et a l. (2002). 
Where only one or two emission lines are present, we 
classify the objects as AGN if broad features or high ion- 
ization lines were present and galaxies otherwise. Note 
that the classifications in most of the optical spectro- 
scopic follow-up wor k in the Chandra Deep Fields, e.g., 
ISzokolv et al. ( 200-1), are based upon both X-ray (lumi- 
nosities and hardness ratios) and optical properties. We 
do not consider the X-ray properties in classifying these 
galaxies as we are using them to construct our X-ray 
prior. 

There are 29 galaxies whose emission-line ratios are 
consistent with starbursts and/or normal galaxies rather 
than AGN. We further exluded two of these sources since 
they were found at off- axis angle in excess of 8' where the 
spectroscopic completeness considerably worsens relative 
to the center of the survey (Szokoly et al. 2004). This 
is thus our "spectroscopic CDF-S sample" with a total 
of 27 galaxies. Note that we do not classify an anal- 
ogous "spectroscopic CDF-N sample" as those optical 
spectra are not available, and also are generally not flux 
calibrated (Ba rger et al.ll2003|h which would prevent line 
ratios from being calculated. 
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Fig. 1. — X-ray hardness ratio HR vs. soft X-ray (0.5 - 2.0 
keV) luminosities for CDF-S sources with spectroscopic redshifts 
and classification. Galaxies are plotted as black squares, AGN1 
are red diamonds and AGN2 are blue circles. All luminosities are 
calculated assuming H = 70 km s" 1 Mpc" 1 , f2 m = 0.3 and fl A = 
0.7 . 



Multi-parameter Classification: {fx/ fa), HR, Lx 

Even with optical classification there is still the 
risk of AGN contamination, particularly for type-2 
AGN where the AGN may be obscured and/or the 
spectral aperture encompasses the e ntire galaxy (e.g., 
iMoran. Filippenko fc Ghornokl 120021) . Therefore, we 
used a second selection approach based on the statistical 
properties of the sources. 

For this statistical analysis we concentrated on the 
CDF-S multiwavelength data since again we had ac- 
cess to the optical spectra of potential counterparts, 
and could therefore most reliably identify control sets 
of galaxies and AGN. We investigated the distributions 
of the X-ray hardness ratios and 1.4 GHz, K-band, R- 
band and soft X-ray luminosities and inferred that a 
good discriminator between AGN and galaxies was based 
on the X-ray/optical luminosity ratio, the X-ray hard- 
ness rati o and X-ray lum inosity as had been found pre- 
viously ijHasinger 1 120031 and references therein). The 
ratios Lx/L\a gh and Lx/L^ also show promise and 
will be explored in future work. 

In Figure ^ w e plot the hardness ratio vs soft band X- 
ray luminosity of all CDF-S sources with good spectro- 
scopic redshifts and classifications. Different symbols are 
used for the various optical spectroscopic classifications 
(type 1 AGN, hereafter AGN1; type 2 AGN; hereafter 
AGN2 and normal galaxies). It is evident that the AGN1 
and galaxies have similar hardness ratios and AGN1 have 
significantly higher X-ray luminosities. AGN2 span a 
larger range in X-ray luminosity than AGN1 and on av- 
erage have harder spectra than AGN1 and galaxies. Both 
of these effects are expected since AGN2 are generally X- 
ray-obscured sources. A brief synopsis of the technique 
follows. 

In most cases the errors on the hardness ratios are large 
for sources in the galaxy regime as there were often only a 
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Fig. 2. — The redshift distribution of the Bayes-selected galaxy 
sample. 



few X-ray counts (i.e., logLx < 42 and HR < 0) and, in 
fact, many sources in this regime only have upper-limits 
for HR. Our selection approach takes into account the 
errors on the hardness ratios and the parent distribution 
properties discussed above. This approach is described in 
Appendix A, and we also restricted our sample to sources 
with redshifts < 1.2 (since the redshift estimates for the 
sources become increasingly uncertain at higher redshifts 
and our two redshift bins should be of comparable size) . 
This resulted in 74 CDF-S and 136 CDF-N sources being 
classified as galaxies. In Figure we show the redshift 
distribution for the (Bayesian) galaxy sample. 

3. LUMINOSITY FUNCTIONS 

The standard method of co nstructing binne d lumi- 
nosity functions discussed in iPage fc Carreral l)200ffl 
was used in this paper to calculate the soft X- 
ray luminosity function for normal galaxies (se e also 
IScnmiddlTflBl IMivau. Hasiuger Schn^l2T)M l20fil 
For each redshift and X-ray luminosity bin, the 
source density, <f>, can be estimated from cf> ~ 

^(JwT SzZf L) ( dVc l dz "> dzdL ^ 1 > where N is thc num " 

ber of sources in the bin, L m i n and L max are the min- 
imum and maximum X-ray luminosities of each bin, 
and dV c /dz is the total comoving volume per redshift 
interval, dz, that the CDF surveys can reach at each 
luminosity L. For a given luminosity Ljf, z m i n {L) is 
the minimum redshift of the bin and z max {L) is the 
highest redshift possible for a source of luminosity L 
for it to remain in the redshift bin. The variance of 
the source densities can also be estimated from 8<p ~ 

SN ULl a ,^z2^ (L \ dV / dz ) dzdL )~ lwheie SN is the lfT 

Poisson error for the number o f sources in ea ch bin 
IKraft. Burrows, fc Nousekl Il99ll lOrfrreTsI Il98^ . Thc 
XLF derived from the Bayesian sample is shown in Fig- 
urc|3and the spectroscopic CDF-S sample XLF is shown 
in Figure 0] 

In the case of the spectroscopic CDF-S sample we also 
made an approximate correction for spectroscopic com- 
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Fi g, 3, — Combined CDF-N and CDF- S XLF based on the Bayesian model testing sample (see text for details) along with the 60/^m LF 
from Takcuchi, Yashikawa, & Ishii (2003) converted to the 0.5-2.0 keV bandpass using the Lgo^m/X-ray correlation (based on the galaxy 
sample given in Ranalli et al. [2003]). The left panel shows the 60/xm LF model based on the full IRAS PSCz sample while right panels 
shows the LF based on the "warm" galaxy sample (with S ioo fim / F^o^im < 2.1). The X-ray data in both plots are identical, and evidently 
the warm galaxy LF more closely matches the XL F at each redshift. The XLF wa s computed with redshifts z < 0.5 and z > 0.5. Also 
plotted (as open circles) is the local XLF derived in Schmi dt. Boiler. &; Vogea 119961) . corrected for the factor of ~ 3 overdensity of galaxies 
locally (see text). Curves are plotted for z=0 and the mean redshifts of the two XLF samples (z = 0.27,2 = 0.79), with the z>0 curves 
being calculated assuming (1+z) 2 7 evolution. The dashed curves show the (minor) effect of including an estimate of X-ray k-corrections 
when converting from FIR to X-ray 0.5-2.0 keV luminosities. The y axis units are log (Mpc - 3 dex -1 ). 



pleteness. We used the mean log Fx / F opt value for galax- 
ies (-1.7 as discussed in the Appedix) to compute the 
range in R magnitude corresponding to each XLF bin, 
and took the fraction of redshifts in that interval (us- 
ing all sources in Szokoly et al. (2004), restricted to a 
maximum off-axis angle of 8') as an estimate of the com- 
pleteness. Often in the higher luminosity bins there were 
not a large enough number of sources for the fraction 
to be computed accurately however we found that the 
completeness when R< 23 was consistently ~ 90%, and 
therefore we fixed the completeness at 90% whenever the 
entire higher R range for that bin was < 23. This pro- 
cedure only significantly impacted the first point of each 
XLF, with completenesses of 43% and 73% for the z<0.5 
and z>0.5 XLFs, respectively. 

3.1. X-ray Luminosity Function Evolution 

Form of the X-ray Luminosity Function 

To describe the evolution of the X-ray luminosity in 
more detail it is useful to use a functional form for the 
XLFs we have obtained above. It is then easy to test 
simple models for evolution of the XLFs such as pure 
luminosity evolution (PLE). The choice of the appro- 
priate functional fit should, if possible, have a solid ba- 
sis in physical understanding and in empirical observa- 
tional correlation. After trying various forms we chose to 
use the lognormal distribution. The infrared luminosity 
function of galaxies (IRLF) is well fit by this distribution. 
In addition, the IRLF is known to be directly pr opor- 
tional to the star- formation rate ljKennicuttJH998|) . We 
established that this form of the distribution is a reason- 
able fit to our XLF data. Then, to give the XLF a more 
physical basis, we predict the XLF from the IRLF based 



on the iRanalli. Comastri fe SettH l)2003f) correlations be- 
tween IR and X-ray fluxes as discussed below. 

IRLF comparison 

We use the 60 iim lu minosi ty function of 
lTakeuchT~Y ashika wa. fe Ishiil l|2003|K which as- 
sumes the same fun ctional form for the IRLF as 
iSaunders et alJ l)1990|) . b ut includes several improve- 
ments . Additionally, Takeuchi. Yashika wa. fc Ishrj 
( 2003) adopt the same cosmology used here, allowing 
for direct comparison. We first have to convert to 
the IRLF to appropriate units to match the XLF. 
The lum inosity function conversions were ca lc ulated 
using (iGeorgantopoulus. Basilakos fc Plionisl 119991: 
lAvni fc Tannenbaumlll986|) : 

(l>x(logLx) = / c/) IR (logL IR )$(logLx\ logL IR )dlogL 
Jo 

(1) 

where $(logLx| log Lj R ) is the probability distribu- 
tion for observing Lx for a given Lj R , which we take to 
be a Gaussian, giving 



IB 



<p(logL x ) 



+ oo i 

(p IR (\ogL IR ) — 



2na 



-(alogL JR + b-\ogL x ) 2 , . 

exp — dlogLm (2) 

We calculated the convolution given in eq. 2 numeri- 
cally with the dispersion of a — 0.3, consistent with the 
dispersion in the soft X-ray / L e0fMm correlation. 

In Figures 01 and 0] we show the 60 um lum inosity 
function of iTakeuchi. Y ashikawa. & Ishii (2003), which 
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Fig. 4. — The CDF-S spectroscopic galaxy sample (see text for details). Curves and symbols are as in Figure |3 with the 60 /im "warm" 
galaxy LF model being used here. The XLF shown here were corrected for completeness, which only impacted the first bin of each XLF 
signficantly (see text). 



assumes the s a me fu nctional form for the IRLF as 
ISaunders et alJ II1990D. In F igure El we plot the 60 /im 
LF from iTakeuchi. Yashikawa. fc Ishiil ((2003) converted 
to the 0.5-2.0 keV bandpass using log L . 5-2.0 koV = 
Leo^m — 3.65 (based on the sample galaxy sample and 
X-ray data as given in Ranalli et al. [2003]). 

The XLF figures also show the effect of including pure 
luminosity evolution of the form (1 + z) 2 - 7 (see §4) for 
the mean redshifts in each rcdsh ift bin (sim i lar to the IR 
evolution 3 ± 1 observed in ISaunders et al.l l)1990j) l. and 
k-correcting the X-ray luminosities derived from the IR, 
i.e., all XLF points are based on luminosities computed 
without k-corrections, equivalent to assuming a flat spec- 
trum with an energy index of ~ 1, while the IRLF curves 
where adjusted as discussed below. We assumed a canon- 
ical SED (shown in Figure El see Ptak et al. 1999) 



with a soft thermal plasma component (kT = 0.7 keV 
and 0.1 solar abundances), representing hot gas possi- 
bly associated with superwinds, and an absorbed power- 
law component representing the X-ray binary emission 
(Nh = 10 22 cm~ 2 and an energy index of 0.8). The 
contribution of the hot gas is softer than the binaries 
but the total effect is to give a flat SED. with SFR. The 
0.5-2.0 keV/2.0-10.0 keV ratio was varied with 0.5-2.0 
keV luminosity as observed in RCS (going from ~ 1.7 
at L x = 10 39 ergs s" 1 to - 0.3 at L x = 10 42 ergs s~\ 
derived from the Lq. 5-2.0 keV and L2-10 keV / FIR corre- 
lations). Therefore the typical case is in fact consistent 
with an effective 0.5-10 keV energy index of ~ 1.0, or a 
vF v slope of 0. 

As a consistency check, we computed the hardness ra- 
tios corresponding to these model spectra. The model 
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Fig. 5. — Spectral energy distribution starburst model used to "k-correct" the models shown in Figures [3] and [3] In the model shown 
here L(0. 5-2.0 keV) = L(2-10 keV), which is the case expected for L(0. 5-2.0 keV) = 10 40 ergs s — 1 (see text) and corresponds to an effective 
energy index of 1.0. The dashed and dot-dashed lines show the contribution of the soft and hard spectral components, respectively. 
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HRs ranged from —0.8 < HR < —0.3 which can be con- 
trasted with the observed HR ^< —0.5 although for the 
vast majority of our sources there was no significant de- 
tection in the hard band. In addition, when the observed 
hard/soft X-ray luminosity ratios are correlated with soft 
X-ray luminosities in RCS, the scatter is very large and 
there is no obvious trend as implied by the correlations 
with FIR. However, since our assumed models are rela- 
tively flat, the k-corrections are minor as can be seen in 

F igure Bl 

Tak euchi. Yashikawa. fc Ishiil l)2003|) note that the 60 
/im LF binned from the full IRAS PSCz sample differs 
from the LFs derived from "warm" and "cool" galaxy 
samples, where "warm" is defined as Sioo/jm/ S&Oum < 
2.1. This flux ratio corresponds to a temperature of 
~ 25K. In Figure El we show both the full and warm 
galaxy sample LFs and evidently the XLF is better 
matched the warm XLF. The XLF is mostly sampling 
L > L* galaxies, and therefore this is due the bright- 
end LF slope being flatter (a = 0.625 as compared to 
a = 0.724 in the log-normal LF parameterization used 
here). This may be due to the fact that the galaxies in 
our sample have relatively high implied star formation 
rates (see below) which results in high dust tempera- 
tures, however we note that AGN activity will also re- 
sult in higher dust temperatures a nd X-ray luminosities 
ijMilev. Neugehanier fc SoTferHTflSl . 

Radio LF Comparison 

The local radio luminosity function is in excel- 
lent agreement with the FIR l uminosity function 
(jCondon. Cotton fe Broderickl 12002). and therefore the 
z=0 FIR model LF is equivalent to the local radio 1.4 Ghz 
LF. This is not surprising considering the tight FIR/1.4 
Ghz correlation. 

Ha Luminosity Function Comparison 

The Ha luminosity is a traditio nal tracer of the m assive 
star-formation rate in galaxies (I Ke nniciittl [19831 . The 
Balmer line strengths are proportional to the number 
of ionizing phot ons from youn g massive stars embedded 
in HII regions ()Zanstralll927|) . Since the Balmer lines 
lie in the red part of the visible spectrum they are not 
as hypersensitive to dust as is the UV continuum. Un- 
like the FIR, Ha can be observed to high-redshift. The 
Ha LF is foun d to evolve by up to an order of magni- 
tude by z = 1 lIGlazehrook et al .lll999l lYan et alJlTflM 
Hopkins, Connolly, fc Szalavl l2000|) . Here we compare 
the XLF to the Ha LF since X-rays also are primarily 
sensitive to the massive star formation rate. 

We converted the older Ha LFs to the now stan- 
dard cosmology adopted for this paper in the follow- 
ing way. Ha luminosity functions were calculated as- 
suming qo — 0.5 and Q\ = 0, with various values 
for Hq. Converting for differences in Hq is straightfor- 
ward. However, the cosmological constant introduces a 
rcdshift-dcpcndence to the luminosity function transla- 
tion. One should return to the original data and re- 
calculate the luminosity function but this requires the 
measured flux values, the redshifts, and detailed knowl- 
edge of the sensitivity which are generally not published 
with the luminosity functions. As a zeroth-order correc- 
tion we evaluate the differences in comoving volume and 



luminosity distance at the median redshift for each lu- 
minosity function, and make the appropriate bulk shifts. 
We may then compare the Ha luminosity functions with 
our X LF. For example, the lHopkins. Connolly, fc Szalavl 
l|2000ft Ha LF has a median redshift of z = 1.3. The 
luminosity distance (D^) for an Qm = 0.3, f^A = 0.7, 
Hq = 70km s _1 Mpc -1 Universe is 46% larger for z = 1.3 
than for the qn = 0.5, Ho = 75 km s" 1 M pc -1 cos- 
mology of Hopkins. ' Connolly, fc Szalavl ((2000). We thus 
shift their Ha LF by a factor of 2.13 to higher lumi- 
nosity. The differential comoving volume (dV c /dz) at 
z = 1.3 is 3.82 times larger for our cosmology, thus the 
iHookins. Connolly, fc Szalavl l)2000Tl Ha LF normaliza- 
tion must be shifted down by this factor. 

However, these shifts will not allow for cosmology- 
dependent changes in the shape of the luminosity func- 
tion (due to the variance in the median redshift of the 
various bins). Since we have all the required data for the 
XLF, we thus also calculate our XLF in a [qo = 0. = 
0,-ffo = 70km s _1 Mpc -1 ] cosmology to ensure that in- 
ferred differences in luminosity functions are not just ar- 
tifacts of our zeroth-order cosmology correction. 

In Figure we show Ha luminosity functions covering 
the full redshift range of interest (z « 0-1.3), for both the 
£!a = 0.7 and fl\ = cosmologies. At low er redshift, the 
compa rison Ha LFs include the work of iGallego et alJ 
(1995) using the Universidad Complutense de Madrid 
(UCM) sur vey for emission-line objec ts (z < 0.045) and 
the work of lTresse. fc Maddox I l)1998t) using the Canada- 
France Redshift Survey (CFRS; median z sa 0.2). At 
higher redshift, Ha shifts into the observed near-infrared; 
the comparison sa mple at z « 0.7 is t he VLT ISAAC- 
observed sam ple oflTresse etaLl l)2002t) . The extinction 
corrections in lTresseetTalJ(|200"2fl are not as reliable as 
the lower-redshift data due to the difficult of measur- 
ing H/3 in the NIR spectra. Finally, at even higher 
redshift, there are two Hubble Spa ce Telescop e Near In- 
frared Camera (NICMOS) studies: lYan et all ()1999fl and 
IHookins, Connolly, fc Szalavl l)2000j) . These two studies 
have median z ss 1.3 and are not extinction-corrected. 

In Figure El the Schec ht er fun ct ion fits to the Ha 
LFs fro m IGallego et al. I (I1995D, ITresse. fc Maddox I 
fiflSa. IHonkins Connolly, fc Szalavl (I200TF a^d 
ITresse et all (I2002T) are shown. Note that the re- 
sults of lYan et all fl 999) are consistent with those 
given in Hopk ins. Connolly, fc Szalavl l)2000?) and 
were included in the Sc h echter function fits in 
IHookins. Connolly, fc Szalavl IpfflM . The Ha lumi- 
nosities are converted to a SFR function by multiplying 
the Ha luminosi ties by 7.9 x 10~ 42 M Q yr _1 ergs -1 s 
l)Kennicuttl I1998D . The XLF based on the combined 
CDF-S + CDF-N Bayes-selected sample shown in Figure 
13 was converted by multiplying the X-ray luminosities 
by 2.2 x 1O~ 4O M yr" 1 ergs" 1 s (RCS). Note that the 
extinction-corrected z ~ 0.7 Ha LF overlaps with the 
uncorrected z ~ 1.3 LF, implying that the extinction 
correction (a factor of ~ 2) is of the same order as the 
evolution in the Ha LF between these redshifts. 

The most striking feature of the XLF / Ha LF com- 
parison is that the two sets of luminosity functions have 
different shapes, and more specifically different bright- 
end slopes. As with the FIR comparison the XLF data 
prefer a flatter bright-end slope than is the case in the Ha 
Schecter function fits. The bright-end of the Schechter 
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Fig. 6. — Comparison of Ha and X-ray luminosity functions (based on the combined CDF-N and CDF-S Bayes-se lected sample as shown 
in Figure[3J after both have been rescale d to star-formation ra tes. Th e Ha points were taken fromlTresse et al.l fiMffi z ~ 0.7), an d Schcchtcr 
functi on fits to Ho LFs we re taken fro m Gallcgo^e^alJ 1199^ z ~ Q.) ITresse. fc Maddox 1 11993, z ~ . 3), a nd Hopki ns. Connolly, fc Szalavl 
120001 z ~ 1.3). Note that lHopkins. Connolly, k, Szalavl 120001 included the Ha data from lYan eTal] il99Gft in their fits. The Ha data are 
listed in the legend in order of increasing redshift; the z < 0.5 data and curves are plotted in blue and the z > 0.5 data and curves are 
plotted in red. The curves are plotted only at the range in luminosity (and hence SFR) over which the fits were performed. The left panel 
shows the X-ray luminosity function calculated for the cosmology H$ = 70 km s — 1 Mpc -1 , Q m = 0.3 and f^A = 0.7, and the published Ha 
values were converted from the cosmologies used (Hq = 50 km s — 1 Mpc -1 , 50 = 0-5 m most cases) by scaling the luminosities by and 
the normalization by dV c /dz~ 1 evaluated at the mean redshift of each survey (see text for details). The right panel shows the XLF rebinncd 
for a cosmology of Ho = 70 km s -1 Mpc -1 and qo = 0.5, in which case the published Ha LFs are simply rescaled by Hq 2 in luminosity 
and Hq in no rmalization. Both cases are shown h ere since the cosmology correction applied in the first case is only an approximation (out 
of necessity). Hopkins, Connollv. & Szalav 12000) gave two Ha LFs based on a "conservative" sample where the Ha emission was securely 
measured and a sample in which less secure Ha line identifications were also considered (i.e., the "true" Ha LF at this redshift is between 
these two extremes). 



function is an exponential cut-off, and therefore the only 
"tuning" that can be done to fit the bright end of the 
LF is adjusting L*. Accordingly the closest agreement 
in shape between the X-ray and HaLFs is with the Hop- 
kins et al. "conservative" Schechter fit, which had the 
largest value of L*. However, note that the shape of the 
Ha LF fits should only be considered to be representa- 
tive since the Schechter function fit parameters are not 
usually well constrained (and in addition are often cor- 
related; see Hopkins et al. 2000). The relevance of the 
form of the luminosity functions is discussed below. 

3.2. Towards a Physical Understanding of the 
Luminosity Functions associated with Massive Star 
Formation Indicators. 

The galaxy luminosity functions for cosmic star forma- 
tion indicators have either: (I) a Schechter distribution 
or (II) a lognormal distribution and we discuss each of 
them in turn below. 

Schechter Distribution 

The Schechter distribution was derived to fit the lumi- 
nosity function for normal galaxies. It is associated with 
the hierarchical build up of luminous baryonic mass, in 
dark matter potential wells, to form the luminous galaxy 
distribution function. It is obvious that if: (1) the initial 
stellar mass function(IMF) is approximately universal; 
(2) the (massive) star formation rate per unit luminosity 



is roughly constant and (3) the Ha is not significantly 
obscured by dust then the HaLF should roughly reflect 
the parent Schechter LF of the galaxies. 



Lognormal Distribution 

Lognormal distributions are associated with physical 
systems that depend on many random multiplicative pro- 
cesses (and many associated random multiplicative prob- 
ability distributions). A good example is the complex 
physics leading to relatively simple and robust proba- 
bility density functi ons in a multi-phase, tu rbulent self- 
gravitating system iWada fc Normanll200l . It is not 
surprising a lognormal distribution results for complex 
star forming systems with, for example, either: (1) repro- 
cessing of the radiation by dust into the IR,to make the 
IRLF or; (2) evolutionary processes of massive stars in 
binaries accreting matter to produce the X-rays of the 
XLF. One way to think of this is to imagine, for example 
in the x-ray binary emission case, all the time ordered 
processes that a piece of matter undergoes on its way to 
emit an X-ray photon from accretion processes onto com- 
pact objects in binary systems. Subsequently, the pho- 
ton may also be reprocessed on its way to the observer. 
To make this argument more explicit let us assume that 
the Schechter function discussed above for the HaLF is 
dependent on the variable X(SHa) given by: 
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X(SHa) = x{star)x(IMF)x(U)x(f B ) = *[[x(i) (3) 

i=i 

where the four x(i) functions on the right are related to 
star formation, the IMF, the cosmological parameters in- 
cluding fi m , and the baryonic fraction respectively. This 
schematic model is meant to indicate that the HaLF de- 
pends on a few variables. 

The physics behind the production of X-ray luminos- 
ity are more complicated. For example for the X-ray 
luminosity variable X(xray) we could write: 

X(xray) = X(SHa) x(binary) x(compact) x(high mass) 

x (accretion) x(metallicity) x(HI)... = YliLi x (i) 

where here the, N, x(i) functions on the right side de- 
pend not only on the four x(i) variables that we have 
included in our schematic representation of the Schectcr 
Ha LF but on a total of, N, x(i) functions describing the 
formation of binaries, compact objects in the binaries, 
high mass companions, physics of accretion, mctallicity 
and HI column ... respectively. Here, the point is that 
the production of X-rays depends on many more vari- 
ables, N » 1, in an approximately multiplicative fash- 
ion. Consequently, the natural distribution for large N 
multiplicative variables is the lognormal. 

Similar ideas apply to the IRLF where, again, we can 
write: 

X(IR) = X(SHa) x(dust) x(transfer) x(metals) 
x(SN:dest/form) x(molcloud) x(m,i) x(dyn)... 

= n£i*(o 

where again we have, in addition to the basic Schecter 
variables, additional x function variables describing the 
physics of dust, radiative transfer, supernovae destruc- 
tion and formation of grains, the molecular environment, 
the triggering of star formation by mergers and interac- 
tions and other dynamical instabilities (such as bars, and 
oval distortions)., respectively. 

In summary, our argument is that the lognormal dis- 
tributions, such as the IRLF, XLF and RadioLF are 
functions of many complicated physical processes, as dis- 
cussed above, and the Schechter type functions, such as 
the HaLF and the K-band LF depend on relatively sim- 
ple and relatively few processes. 

Future Work 

Detailed models could be developed using the ideas 
discussed above but with more attention to the detailed 
physical processes. Possibly, the lognormal distributions 
may be associated more with starbursting normal galax- 
ies rather than normal quiescent galaxies. It is likely 
that starbursting events are triggered in some way, by 
external galaxy merging and interactions or by internal 
dynamical instabilities such as bar formation. The ran- 
dom nature of such triggering processes may naturally 
lead to the type of multiplicative probability chain that 
would produce lognormal distributions. The more nor- 
mal star formation mode may be reflected by the Ha. It 
would then be very interesting that normal modes and 
the triggered modes contribute very roughly the same 
amount of the cosmic SFR. For a detailed discussion of 
the merits of Ha v ersus IR star formation rate indicators 
see lKewlev eTafl l|200^ . 



4. COSMIC STAR FORMATION HISTORY 

In Figur e [71 we show the cos mic star formation history 
taken from lTresse et all l|2002f) . with X-ray points added 
at z=0.27 and z=0.79 based on the combined CDF-N+S 
Bayes sample. We chose the Bayes sample since the 
XLF is in better agreement with the FIR prediction 
(see below). We computed the X-ray luminosity den- 
sity by integrating Q(logL)LdlogL and rescaling by 2.2 x 
10- 40 M o yr" 1 erg s" 1 (RCS). This resulted in SFR den- 
sity estimates of 2.0 x 10~ 2 and 6.0 x 10~ 2 M yr^Mpc" 3 
at two redshifts. 

Note that the integration of an LF tends to be domi- 
nated by the contribution of the LF near L* . Here we pri- 
marily only sample galaxies brighter than L* , and there- 
fore this estimate is a lower limit to the true SFR. How- 
ever we also integrated the rescaled FIR LF (the curves 
in Figure 03 including k-corrections) at those redshifts, 
resulting in SFR = 0.022 and 0.055 M yr^Mpc" 3 . 
The two approaches resulted in consistent SFR estimates 
however the errors are large on the XLF data points, 
particularly on the lowest luminosity points which dom- 
inate the luminosity density sum. Therefore, we esti- 
mated the error on the SFR by summing the upper and 
lower error bounds on the XLF data points. This re- 
sulted in errors in log SFR of ~ 0.16 and ~ 0.31 at 
z=0.27 and 0.79, re spectively. In this fi gure the SFR den- 
sity estimates from lTresse et all l|2002f) were corrected to 
h = 0.7, fix = 0.7, fl M = 0.3. 

The X-ray SFR points are consistent with an evolution 
of the SFR for < z < 1 of SFR oc (1 + z) 2J (i.e., an 
increase by a factor of ~ 2.5 in SFR density between z 
~ 0.3 and z ~ 0.8). The Tresse et al. SFR history plot 
includes UV and extinction-corrected Ha values. The 
z ~ 0.3 XLF SFR estimate is more consistent with the 
Ha SFR estimate and the z ~ 0.8 XLF SFR estimate is 
intermediate between the Ha and UV SFR values. How- 
ever since the error on the z ~ 0.8 XLF SFR is large 
it is consistent with either the UV or H a SFR. We also 
plotte d the z ~ SFR predicted from the lSaunders et all 
(1990) 60 /im LF (solid line) with an evolution of (1+z) 3 , 
as well as th e 60 ^m SFR predicted by integrat ing the 
60 fj,m LF in lTakeuchi. Yashikawa. & Ishiil l|2003L dotted 
line). In both cases the 60 /mi/SFR conversion of 2.6 x 
10" 10 M Q yr-^Leo/Lo]- 1 (whe re L m /Lp, is the 60 /mi 
lumin osity in solar units) from | R ^w^nT{ obinson^^^ 
(1997) was used. The iTakeuchi. Yashikawa. fc Ishil 
( 2003) 60 /im LF has a somewhat steeper faint-end slope, 
which resulted in a factor of ~ 2 larger luminosity den- 
sity. 



5. DISCUSSION 

Our studies have focused on establishing the X-ray de- 
rived cosmic star formation rate and the normal star- 
forming galaxy X-ray Luminosity functions up to red- 
shifts of order unity. Mild evolution is seen. The results 
presented here open up a new waveband for such star for- 
mation studies with the significant advantage that obscu- 
ration effects are not present. The Bayesian techniques 
developed in this paper for the galaxy selection have a 
wide range of applicability for further such studies. We 
now discuss important aspects of our results in more de- 
tail below. 
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Fig. 7. — Compilation of SFR densities fromfTrcsse il al. (2002), including the X-ray points for z=0.27 and z=0,79 from this work. The 
X-ray points are shown with black filled squares. The triang les represent the Ha SFR values from lGallego et alj 119951) at z ~ O.IGronwalll 
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11999) at z~ 0.05, and IHopkins. Connolly, fc S zalav (2000) at z~ 1.3. The empty s quare represents the UV-selected Ha IPascual et aT, 
d200ll) value at z~ 0.24. The filled circles show the SFR densities from lTresse et al.l 120021) . The star gives the UV-selected z=0.15 SFR 
density from Sullivan ct al. (2000). The open circles give 2800A CFRS points from L illy et alj 11993) . without dust correction. We also 



plot the SFR history based on the 60 ^im LF luminosity densi ty at z ~ 0, incl uding 
computed from Saunders et al. (1990, solid line) and ffi^euchj/Yashikawa 



evolution. The 60 ^m luminosity density was 
dotted line). 



5.1. Normal Star-Forming Galaxy X-ray Luminosity 
Functions 

We have derived the first X-ray galaxy luminosity func- 
tions at the redshifts of ~ 0.3 and ~ 0.7. The high- 
luminosity bins for both the z<0.5 and z>0.5 Bayesian 
sample XLFs tend to be flatter than the predictions 
based on rescaling the FIR LF, although there is bet- 
ter agreement when the XLF is compared to the "warm" 
galaxy FIR LF rather than the full PSCz sample FIR 
LF. This implies that there is some AGN contamina- 
tion in the sample. The spectroscopic CDF-S sample 
XLF also appears to be somewhat flat, which may indi- 
cate that absorbed and/or low-luminosity AGN are being 
missed in the optical spectra, although the small num- 



ber of galaxies in that sample results in large statistical 
errors. Near IR spectra may be more revealing (e.g., 
IVeilleux. Sanders fc Kin3 ll999L This overall agreement 
with the FIR prediction suggests that X-rays may be a 
useful tool for discriminating AGN2 and galaxies in sur- 
veys, due to the ability of X-rays to penetrate columns 
of 10 22 ~ 23 cm" 2 as typically observed in Compton-thin 
AGN2. Finally, \ogLx < 40.5 bin from the Bayesian z 
< 0.5 XLF is most consistent with the local (full sample) 
FIR L F and the z=0 XLF from lSchmidt. Boiler. & Vosei 
(1996), implying no evolution. However the mean red- 
shift of galaxies with logLx < 40.5 is only 0.16, implying 
evolution only on the order of ~ 50% in luminosity. 

5.2. Star Formation Rates from X-ray Studies 
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Detailed physical modeling of the X-ray emission 
from galaxies and its application in measuring cosmic 
star formation history has been co nsidered recently by 
iCavalier ^ Giacconi fc M encil i2000(l for the hot ga s com- 
ponent. iGhosh fc White! l|200lD and iPtak et alJ (pCIOlh 
for the X-ray binary component, and for hi gh mass bina- 
ries by iGrimm. Gilfanov fc Sunvaevl l)2003[) . The spec- 
tral energy distri bution (SEP) for starburs t galaxies has 
been modeled bv lPersic fc Rephaelil l|2002j) . The poten- 
tially dominant contribution from Ultra-Luminous X-ray 
po int sources (ULXs) ha s been analyzed in a recent paper 
hv lColbert et {200% 

The physical reason for the correlation of both hard 
and soft X-ray luminosity with SFR is worth studying 
because it seems to constrain the overall physics of mas- 
sive star formation and evolution. One possibility is that 
the (accretion-powered) core-collapse supernovae provide 
the hot 1 keV component and the accretion-powered 
high-mass X- ray binaries dominating the point source 
contribution (Grimm, Gilfa nov fc Sunvaevll2"00l . How- 
ever, the detailed physics of how the relative branch- 
ing ratios to SNe and to HMXBs for the evolution of 
the massive star population can conspire to produce the 
flat SED is not yet clear. We argue below that not 
only must the single star IMF be fixed but the bivariate 
binary star IMF(M 1 ,M 2 ) must also be approximately 
constant. For a flat SED, what is required is that the 
massive stars ending in Type II SNe powering the 1 
keV hot gas component of the ISM and the HMXBs 
for the evolution of massive stars in binaries must be 
giving approximately equal energy per octave. Detailed 
studies of nearby galaxies and detailed population syn- 
thesis models help significantly ij Zezai^fc Fabbiano|[2002j: 
Georgantopoulos. Zensas fc Ward! 120031 ISwartzetltfl 
2003HColbert et al.ll2003HSipior et al.ll2003j) . In general, 
the X-ray luminosity of a galaxy can be related to the 
mass of the galaxy, M and its star formation rate, SFR 
by: 

L x =VmM g + V*M* (4) 

where the 77's are linear efficiency factors. The first 
term is due to the old low mass X-ray binary popula- 
tion(LMXBs) and the second term is due to the high 
mass binaries (HMXBs) and the hot gas powered by 
SNe. In the Chandra band (0.3-8 keV), the X-ray point- 
source luminosity is proportional to the mass and SFR 
with coefficients t)m = 1-3 X 10 29 ergs s" 1 M^ 1 and 
p* = 0.7 x 10 39 ergs s" 1 (M@ yr" 1 )" 1 llColber't et all 
2003). In the soft band (0.5-2.0 keV) the star forma- 
tion dominates for galaxies with Lx > 10 39 ergs s _1 (as 
indicated in X-ray observ ations of the composite emis- 
sion of starburst g alaxies; iDahlem. Weaver fc Heckmanl 
(1998); IPtak et alJ IH999DV It was, therefo re, not sur- 
prising that lRanalli. Comastri fc Settil l|2003|) found that 
the soft X-ray luminosity of starburst galaxies is propor- 
tional to the SFR alone, specifically with t\m ~ and 
?7* = 4.5 x 10 39 ergs s" 1 (M yr -1 ) -1 . 

The extinction-corrected Ha star formation history 
and the X-ray star formation history formally agree 
within the errors, although the X-ray SFR values are 
lower than the Ha SFR values. This is probably due in 
part to the different sample selection and the different 
corrections applied to calculate the SFR. The revision 



to the 60 /im LF given in lTakeuchi. Yashikawa. fc Ishiil 
(2003) resulted in a factor of ~ 2 increase in the 
60 luminosity density relative to the value given in 
iSaunders et alJ <)1990j) . It is likely that future work on 
both the Ha LF and XLF will similarly result in adjust- 
ments to the SFR estimates. 

The recent comprehensive study of lCohenI l)2003[) using 
spectroscopy of the [OII]3727 A line emitted by galaxies 
found in the 2 Ms exposure of the CDF-N shows that the 
Ljf-SFR rate is similar locally and at redshifts of order 
unity, which is consistent with what we assume here. 

The basic ingredients of a workable model are fairly 
clear: 

(1) Constant IMF(Mi,M 2 ). 

A constant two-parameter IMF(B(Mi, M 2 ),S(Mi))for 
binary stars(B) and single stars(S) as a function of M\ 
and M 2 . T his is probably an important im plication of 
the work of Ranall i~Comastri fc Settil l)2003|) and if the 
general consistency here between the x-ray star forma- 
tion results and the star formation indicators. 

(2) Normal Stellar Evolution and Binary Star Evolu- 
tion. 

Straightforward stellar evolution and binary evolution 
to produce roughly constant branching ratios for high- 
mass X-ray binaries (HMXB) and core collapse SNell. 
This is probably satisfied since we expe ct no surprises in 
stella r evolution and binary evolution l)Ghosh fc White! 
12001 . 

(3) Consistent ISM properties. 

It is also necessary to have at least very roughly con- 
stant properties for the interstellar medium(ISM). Over 
the whole galaxy ensemble, no large variations in the 
properties of the ISM including: dust-to-gas ratio, grain 
properties, metallicity, dense molecular environments 
and galactic winds. However, the dispersion of ISM prop- 
erties from galaxy-to-galaxy is known to be large. For ex- 
ample, there are large well-known variations in the grain 
properties in the LMC, the SMC and the Galaxy. 

(4) Consistent Dynamical Properties. 

The population of normal starbursting galaxies we are 
studying can have additional dynamical and environmen- 
tal constraints due to the amount of merging, interac- 
tions, barred driven inflow... 

The first two points (1,2) give a way to understand the 
rough correlations and the last two points (3,4) give the 
dispersion and the overall lognormal shape to the radio, 
IR and X-ray LFs. 

The star formation history derived from the X- 
ray data using the empirically derived relation be- 
tween X-ray luminosity and sta r formation rate from 
iRanalli. Comastri fc Sett] i|2003ft is in rough agre ement 
with other methods although, as discussed by iHogd 
l)2004|) there is a large dispersion in the results. 

5.3. Hidden Active Galactic Nuclei 

An important issue is whether the X-ray emission from 
our galaxies could be due to central AGNs hidden in faint 
Type 2 Seyfert nuclei lurking in the ce nter of our galaxies 
iMoran. Filinpenko fc Chornokl l|2002ft . This is not likely 
because Seyfert 2s have a hard spectrum with hardness 
ratio HR > 0, whereas the galaxies in our sample have 
softer spectra. 

A key advantage to using soft X-ray luminosities 
to estimate the X-ray galaxy luminosity function and 
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the SFR is that the soft band tends to be domi- 
nated by starburst emission regard less of whether an 
obscured AGN is present (see, e.g.. iTurner et alJll997t 
ILevenson. Weaver. &: HeckmarJl200lj) . This assists us in 
discerning galaxies from AGN2. 

5.4. Implications for Future X-ray Missions 

Our cosmological star formation studies, extending 
up to z rs 1, require the combination of 1-2 Msec 
Chandra exposures and high-quality optical spectro- 
scopic coverage from large-aperture ground-based tele- 
scopes such as VLT and Keck. There is, however, 
a wealth of information on galaxy evolution available 
through current ground-based wide-field optical surveys 
(e.g., 2dF and SDSS) that is still relatively underuti- 
lized. The small fields of view of Chandra and XMM- 
Newton make X-ray studies of thes e wide-field galaxy 
surve ys quiet difficult (but see e.g., iGeorgakakis et all 
120031 for an example of one high-quality statistical study 
of 2dF galaxies with X MM- Newton). Deeper wide- 
field X-ray surveys would also be complementary to the 
deep wide-field ultraviolet imaging surveys of GALEX. 
The GALEX surveys, currently underway, are prob- 
ing star formation from z — — 2, overlapping with 
the work we present here. Possible future wide-field 
X-ray telescop es include the Wide Field X-ra y Tele- 
scope (WFXT) iBurrows. Burg fc Giaccon il il992t) a nd its 
successors WAXS/WFXT llOhmcarim et al.ll199Sh and 
DUET (Jahodaet al. 2003). 

X-ray studies of galaxies may be extended to truly ex- 
treme redshifts when high-sensitivit y X-ray telescopes 
such as XEUS l)Parmar et alJ HQ03) commence opera- 
tion. XEUS may achieve sensitivities as faint as 4 x 
10~ 18 erg cm -1 and at this level galaxies will dominate 
the X-ray source number counts. This sensitivity will 
allow for the detection of star-forming galaxies to great 
distances (z > 5) and will allow for the probing of X-ray 
emission from star formation to the very epoch of galaxy 
formation in the early Universe. 

6. CONCLUSIONS 

There are five main conclusions of this work. 

1. We have established the X-ray luminosity functions 
of galaxies in the Chandra deep Fields North and South. 
These XLFs have a lognormal distribution. The XLFs 
are completely consistent with both the infra-red lumi- 



nosity functions (IRLFs) and the GHz radio luminosity 
functions (RLFs), which both have lognormal distribu- 
tions. 

2. The evolution of the XLFs is consistent with a pure 
luminosity evolution(PLE) of ~ (1 + z) 2 7 .The more ro- 
bust integral of the XLFs can be transformed to an x- 
ray derived cosmic star formation history. This cosmic 
star formation rate (SFR) is consistent with the (mild) 
evolution in SFR oc (1 + z) 2 - 7 . In general, the over- 
all SFR estimates in the range of redshift from unity to 
the present epoch show a wide scatter although there 
is a general trend of mild evolution consistent with the 
above. To reduce the scatter and better understand the 
evolution of the SFR in this redshift range one needs 
in the future:(l) more sensitive wide area studies such 
as the GALEX mission (2) extensive multi-band studies; 
the XLFs and X-ray derived SFRs now bring in a new 
waveband here and (3) deeper physical understanding of 
the star formation indicators. 

3. The HaLFs have the form of the Schechter luminos- 
ity function which fits the galaxy luminosity functions in 
J and K. The H a distribution is quite different from the 
the lognormal distributions for the XLF, IRLF and RLF 
discussed above. A detailed physical understanding is 
not yet available. 

4. The X-ray derived star formation history gives an 
interesting constraint on the IMF. The X-ray emission 
in the Chandra band is most likely due to binary stars. 
The overall consistency and correlations between single 
star effects and binary star effects in different wavebands 
indicate that the bivariate IMF(Mi, M2) must be univer- 
sally constant, at least at the high mass end, since the 
X-ray observations may be measuring directly the binary 
star formation history of the Universe. 
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APPENDIX 

BAYESIAN STATISTICAL ANALYSIS 

For many of the sources in our sample both the X-ray and optical information is limited. Classifying a source 
based simply on "sharp" cuts in parameter space, such as log Lx < 42 ergs s _1 and X-ray hardness HR < —0.2, 
does not take into account the errors in the measured quantities and the fact that the parent distribution (i.e., real 
galaxies in the Universe) are not so sharply delineated. Here we remedy this situation by using a Bayesian statistical 
approach. Specifically, we assume parent distribution models for measured parameters, 9, and then compute the 
conditional probability, Pm, of measuring the observed values of the parameters (D) by integrating over the parent 
distributions for a given model M. Specifically, 9 = {HR, \ogLx, log Fx/F opt }, where HR is the X-ray hardness 
ratio, Fx is the 0.5-2.0 keV luminosity, and logFx/F opt is the X-ra y /optical flux ratio, log F x /F mlt is given by 
log 10 (i<o. 5-2.0 k R y — R/2.5 — Rq, wh ere Rq = 5.5 for CDF-N sources l)Hornschemeier et all 12003) and Rq = 5.7 for 
CDF-S sources l|Szokolv et al.ll200l . P M is given by P M = J p M {0\D)d0, with 

p M (9\D)=p M (9)p(D\9)/p(D) (Al) 

p(D\6) I p(D) is the probability of observing a given set of parameter values D when the expected (i.e., parent) values 
are 9, which is simply the likelihood function for the data. In this analysis we assume that the errors on each parameter 
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Table Al. Statistical Parameters for the Galaxy, AGN1 and AGN2 
CDF-S Samples 



















Source Type 


Number 


logLx 


°"log L x 


HR 




\ogFx/F opt 


a lo S F x /F apt 


Galaxy 


8/29 


41.0 


0.8 


-0.26 


0.19 


-1.7 


0.6 


AGN1 


25/27 


43.4 


0.5 


-0.48 


0.08 


0.0 


0.4 


AGN2 


26/28 


42.1 


0.9 


0.21 


0.23 


-0.7 


0.6 



Note. — The Number column lists the number of sources of the given type with a detection in the hard band / the total number of 
sources of that type. Sources with upper-limits on the hardness ratio (i.e., without hard band detections) were excluded when computing 
the mean and standard deviation of HR. 



Table A2. Bayesian Classification Results 



Source Type 


Number 


Galaxy 


AGN1 


AGN2 


Galaxy 


29 


25 


2 


2 


AGN1 


27 


1 


25 


1 


AGN2 


28 


3 


1 


24 



Note. — Source Type gives the spectroscopic classification of the sources. 



are reasonably approximately Gaussian, and defer the treatment of the more general likelihood functions for future 
work, pm is the prior, or parent population probability distribution for a sources of type M. We further assume that the 
parent distributions of the parameters are also reasonably approximated by a Gaussian, with the mean and standard 
deviations of the Gaussian set to the observed values for the parameter (see below), pu should be normalized by 
the relative number of sources of a given type that would be detected in the 0.5-2.0 keV bandpass. However, since 
these numbers are not known precisely and the relative number of sources used in the statistical analysis below are 
approximately equal, we set this term to unity here. Therefore Pm is given by 

P TT [ M 1 1 ( g »-^) 2 1 1 r - ft) 2 ! , A0 , 

Pm ^1 m ^ eM ~—^r ] ^K^ l eM —^r ] (A2) 

Here 9i and <7j are the parent distribution mean and standard deviation for the ith parameter, while Xi and Axj are 
the measured value and erro r for that parameter. For more discussion of Bayesian methods used in model testing, see 
Hob son fc McLachlanl l)2003l) and references therein. 

Applying these methods to our sample, we have chosen Gaussian parent distribution functions whose mean and 
standard devia tions (a) are derive d from the subset of X-ray sources which have high-quality optical spectroscopic 
identifications l)Szokolv et al]l2"?)04l) 14 . We further require that the sources have constrained hardness ratios (relevant 
mostly for the galaxies, which were X-ray faint) when computing the mean and standard deviation of the hardness 
ratios. The h istograms of the observed parameters are shown in Figures lA8IA10l The Gaussian parameters are shown 
in Table lATl and the corresponding Gaussian curve is also plotted in Figures lA8l A 101 

Using these assumed models, the observed errors on HR, and an assumed error of 0.25 on logLx and logFx / F opt , 
we computed the probability of each source being a galaxy, an AGN1 and an AGN2. Note that the e rrors on HR 
were also taken to be Gaussiam using the larger of the upper and lower error estimates (derived using lLvonslll99"ll 
iKraft. Burrows, fc Nouse"kl 119911 ; see discussion in Alexander et al. 2003). A source is clas sified as a given type 
when the computed probability for that type is the greatest. As a check, we summarize in Table lA"2l the results for the 
sources with high-quality classifications (i.e., those used to produce the parent distribution models as discussed above). 
Of course in general a control sample is necessary to properly assess the performance of a classification procedure, 
however our sample is too small to allow for this approach. However in future work this test will be performed with 
other data sets (most notably using GOODS data). As would be expected, nearly all of the false galaxy identifications 
were AGN2 and vice- versa. 

GALAXY SAMPLE 

In Table lB3l we give the identification numbers (as given in Giacconi et al. 2002 and Alexander et al. 2003), redshift 
and logLjf values for the galaxies in the Bayesian-selected sample. The sources that are also in the spectroscopic 
CDF-S sample are marked. There are also 3 additional sources in the spectroscopic CDF-S sample that were not 
selected as galaxies by the Bayesian analysis, which have IDs 49, 242, and 519, redshifts 0.53, 1.03, and 1.03, and 
\ogLx values 42.3, 42.3, and 41.9. As might be expected, these were the spectroscopic CDF-S sample galaxies with 
the highest Lx values. 

14 Ambiguous cases with Lx > 10 42 ' 5 ergs s — 1 where classified as AGN1 galaxies 
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Table B3. XLF galaxy sample 



ID # a z L x (0.5-2 keV) 



CDF-S sample: 



1 


0.347 


41.73 


12* 


0.251 


41.65 


40 


0.545 


42.16 


77 


0.622 


42.26 


84 


0.103 


40.38 


95 


0.076 


39.95 


96 


0.274 


40.78 


97 


0.180 


41.16 


98 


0.279 


41.02 


103 


0.215 


40.92 


116 


0.076 


39.87 


170 


0.664 


41.55 


173* 


0.524 


40.92 


175* 


0.522 


41.12 


185 


0.925 


41.51 


186 


1.158 


41.81 


211 


0.679 


41.81 


218 


0.497 


41.29 


224* 


0.738 


41.77 


229* 


0.103 


39.78 


233* 


0.577 


41.01 


236 


0.731 


41.68 


246 


0.690 


41.98 


247 


0.040 


38.64 


249 


0.964 


41.85 


504 


0.541 


41.53 


509 


0.556 


41.59 


511 


0.767 


41.52 


512* 


0.668 


41.52 


514 


0.103 


39.63 


516* 


0.665 


41.47 


521* 


0.131 


39.94 


525 


0.229 


40.49 


534 


0.676 


41.03 


535* 


0.575 


41.43 


536 


0.444 


41.10 


538 


0.310 


40.50 


552 


0.673 


41.63 


553 


0.366 


40.97 


554 


0.225 


40.67 


556 


0.635 


41.33 


557 


0.500 


40.85 


558 


0.585 


41.49 


559 


0.114 


39.78 


560* 


0.669 


41.36 


565* 


0.363 


40.42 


566 


0.734 


41.74 


567* 


0.456 


40.81 


573* 


0.414 


40.61 


575* 


0.340 


40.45 


577* 


0.547 


41.23 


578 


0.969 


41.62 


580* 


0.664 


41.13 


581 


0.799 


41.44 


582* 


0.241 


40.15 


586* 


0.580 


41.05 


587* 


0.246 


40.07 


590 


0.280 


40.36 


592 


1.150 


41.78 


594 


0.733 


41.91 


617 


0.588 


41.54 


619 


0.050 


38.91 


620* 


0.648 


41.15 


621 


0.290 


40.30 


624* 


0.668 


41.17 


625 


1.151 


41.76 


627* 


0.248 


40.12 


628 


0.200 


39.98 


629 


0.410 


40.71 


644 


0.119 


40.08 


646 


0.438 


40.79 


650 


0.223 


40.67 


651 


0.182 


40.57 


652* 


0.077 


39.24 



a ID number for CDF-S sample is the XID as listed in Giacconi ct al. (2002). 
These sources are also in the spectroscopic CDF-S sample. 
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Table B4. XLF galaxy sample 



ID # a z L x (0.5-2 keV) ID # a z L x (0.5-2 keV) 



CDF-N sample: 



3 


0.138 


40.29 


6 


0.135 


40.83 


22 


0.317 


39.89 


46 


0.207 


40.43 


49 


0.296 


41.27 


55 


0.637 


41.30 


56 


0.108 


39.77 


57 


0.375 


40.87 


60 


0.529 


40.78 


62 


0.087 


38.70 


66 


0.333 


41.24 


67 


0.639 


41.39 


69 


0.520 


41.11 


78 


0.747 


41.66 


81 


0.409 


40.77 


87 


0.136 


39.76 


93 


0.275 


42.08 


101 


0.454 


40.81 


103 


0.969 


42.13 


105 


0.319 


39.99 


111 


0.515 


40.64 


114 


0.534 


40.63 


119 


0.473 


40.83 


120 


0.694 


41.27 


121 


0.520 


41.33 


126 


0.779 


40.91 


131 


0.631 


40.62 


132 


0.647 


41.07 


136 


0.472 


40.65 


138 


0.483 


40.79 


141 


0.746 


41.87 


148 


1.130 


41.96 


166 


0.455 


40.50 


169 


0.845 


41.53 


175 


1.014 


41.99 


177 


1.016 


42.14 


180 


0.456 


41.24 


189 


0.410 


41.01 


192 


0.680 


40.80 


197 


0.081 


38.92 


200 


0.971 


41.31 


203 


1.143 


41.27 


207 


0.300 


41.04 


209 


0.510 


41.36 


210 


0.848 


41.50 


211 


0.846 


41.35 


212 


0.943 


41.77 


214 


1.144 


41.37 


215 


1.006 


41.17 


218 


0.089 


39.31 


219 


0.845 


40.96 


224 


0.700 


41.92 


225 


0.290 


40.30 


227 


0.556 


40.64 


230 


1.011 


41.62 


234 


0.454 


40.85 


241 


0.851 


41.83 


244 


0.971 


41.36 


245 


0.321 


39.90 


249 


0.475 


41.30 


251 


0.139 


39.59 


255 


0.113 


40.31 


256 


0.593 


41.09 


257 


0.089 


38.92 


258 


0.752 


40.96 


260 


0.475 


40.72 


264 


0.319 


40.26 


265 


0.410 


40.74 


269 


0.358 


40.36 


274 


0.321 


40.92 


279 


0.890 


41.07 


280 


0.850 


41.12 


282 


0.202 


39.87 


285 


0.288 


41.03 


288 


0.792 


41.49 


291 


0.517 


40.59 


294 


0.474 


41.42 


295 


0.845 


41.31 


296 


0.663 


41.11 


300 


0.137 


39.54 


305 


0.299 


40.15 


308 


0.515 


40.53 


310 


0.761 


40.91 


311 


0.914 


41.33 


313 


0.800 


41.30 


316 


0.231 


40.32 


320 


0.956 


41.38 


326 


0.359 


40.30 


327 


0.913 


41.16 


332 


0.561 


40.70 


333 


0.377 


41.62 


337 


0.902 


41.04 


339 


0.253 


39.95 


346 


1.017 


40.98 


351 


0.940 


41.59 


353 


0.422 


40.80 


354 


0.568 


41.05 


355 


1.027 


41.77 


356 


0.956 


42.20 


358 


0.907 


41.86 


359 


0.902 


41.62 


378 


1.084 


41.40 


383 


0.105 


39.18 


387 


1.081 


41.70 


388 


0.559 


41.05 


389 


0.557 


40.78 


392 


0.411 


40.51 


395 


0.411 


40.77 


401 


0.935 


41.51 


4U4 


0.104 


39.26 


407 


1.200 


41.94 


410 


0.113 


39.30 


414 


0.800 


41.46 


415 


0.116 


39.60 


418 


0.279 


40.35 


425 


0.214 


40.11 


426 


0.159 


39.75 


428 


0.298 


40.47 


433 


1.022 


41.64 


435 


0.201 


39.99 


436 


0.189 


40.38 


438 


0.220 


40.46 


443 


0.231 


40.69 


446 


0.410 


40.56 


450 


0.935 


41.74 


453 


0.838 


42.11 


454 


0.458 


41.48 


458 


0.069 


39.55 


460 


1.084 


42.17 


462 


0.511 


40.44 


466 


0.440 


40.81 


469 


0.188 


40.09 


471 


1.170 


41.90 


476 


0.475 


40.95 


480 


0.456 


41.19 


489 


1.024 


42.24 



a For the CDF-N sample, it is the ID number listed in Alexander et al. (2003). 
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Fig. A8. — The observed parameter distributions for AGN1 sources in the CDF-S. The y-axes are numbers of sources. The curves show 
the Gaussian resulting from taking the mean and standard deviation of the distribution. 
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